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The electronic structure of the molecular compound (TTM-TTP)i3, which exhibits a peculiar intra- 
molecular charge ordering, has been studied using multi-configuration ab initio calculations. First we 
derive an effective Hubbard-type model based on the molecular orbitals (MOs) of TTM-TTP; we set 
up a two-orbital Hamiltonian for the two MOs near the Fermi energy and determine its full parameters: 
the transfer integrals, the Coulomb and exchange interactions. The tight-binding band structure obtained 
from these transfer integrals is consistent with the result of the direct band calculation based on density 
functional theory. Then, by decomposing the frontier MOs into two parts, i.e., fragments, we find that 
the stacked TTM-TTP molecules can be described by a two-leg ladder model, while the inter-fragment 
Coulomb energies are scaled to the inverse of their distances. This result indicates that the fragment picture 
that we proposed earlier [M.-L. Bonnet et al: J. Chem. Phys. 132 (2010) 214705] successfully describes 
the low-energy properties of this compound. 

KEYWORDS: multi-orbital molecular compound, intra-molecular charge ordering, ab initio calculation, 
band calculation, fragment decomposition 



A rich variety of physical properties have been explored 
in molecular solids since the realization of the first molec- 
ular metal, TTF-TCNQ, and in the last decade, charge- 
ordering (CO) phenomena and related issues have attracted 
much attention. 1 ' 2 ' Despite the complexity of the constituent 
molecules, a simple picture works well for many compounds; 
the low-energy electronic properties are determined by a sin- 
gle frontier molecular orbital (MO). Therefore each molecule 
can be regarded as a single site, and then the tight-binding 
and Hubbard-type models based on this picture successfully 
describe their physical properties. 2 ' 3 ' However, it has recently 
been recognized that such a simple single-MO approxima- 
tion is not sufficient for describing the electronic structures of 
single-component molecular metals M(tmdt)2 [Al — Ni, Au, 
etc.]; 4-6 ' there are cases where multi-orbital effects reflecting 
the respective MOs induce unexpected novel physical proper- 
ties. 

In the present paper, we focus on a quasi-one-dimensional 
molecular crystal (TTM-TTP)l3. 7 ' 8 ' Since the formal charge 
of the TTM-TTP molecule is +1 owing to the monovalent 
counterion 1^, this compound was considered to be a half- 
filled system and the possibility of a genuine one-dimensional 
Mott insulator with paramagnetic spin properties was dis- 
cussed. 8 ~ 10 ' However, it has also been pointed out that this 
system exhibits a transition toward a nonmagnetic state at 
around 80 K. 1 '~ 13 ' On the basis of detailed experimental anal- 
yses of this low-temperature phase by Raman scattering 14 ' 
and X-ray measurements, 15 ' a new type of CO state has been 
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proposed. This state is called the "intra-molecular CO" (ICO) 
state, where the inversion center on the middle point of the 
TTM-TTP molecule is lost and the charge is disproportion- 
ated within each molecule. This state cannot be described by 
the simple single-MO picture by its nature. 

A key to understanding the electronic state of this com- 
pound, as shown in our previous theoretical work, 16 ' is that the 
singly-occupied-molecular-orbital (SOMO) and the second- 
highest-occupied-molecular-orbital (HOMO-1) of the ionic 
[TTM-TTP] + molecule are close in energy. The chemical rea- 
son of this quasi-degeneracy can be understood in terms of a 
three-fragment model. 16 ' The SOMO and HOMO-1 are basi- 
cally described by two of the fragments, as in the following. 
Since the TTM-TTP molecule itself has an inversion center, 
the resulting MOs can be classified into gerade (g) and unger- 
ade (u) MOs. The SOMO (u MO) and HOMO-1 (g MO) 
are shown in Fig. 1, where the SOMO exhibits a bonding 
character between the left and right fragment MOs while the 
HOMO-1 has an antibonding character, i.e., 



(-<PL + <Pk) 



(1) 



where </?l and </?r are the left (L) and right (R) fragment 
MOs. 16 ' The purpose of the present paper is to construct a 
minimal model that can describe the ICO state. The con- 
tribution of the center (C) fragment MO can be neglected 
for this purpose, since the C fragment is deep in energy. 16 ' 
We construct a two-orbital Hubbard-type model based on the 
SOMO and HOMO-1. The magnitudes of the intra-molecular 
and inter-molecular interactions are estimated from ab ini- 
tio calculations with a complete-active-space configuration- 
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Fig. 1. (Color online) SOMO and HOMO-1 of the ionic [TTM-TTP]+ 
molecule obtained by the restricted open-shell Hartree-Fock procedure, 16 ' 
and a schematic energy diagram based on the left and right fragment MOs. 



interaction (CAS-CI) method. 17 ' By transforming the model 
into the fragment picture, we clarify that the stacked TTM- 
TTP molecules can be described by a two-leg ladder model. 

The ab initio calculations were performed using the Mol- 
cas7 package 17 ' with the basis set: S(7s6pld)/[4s3pld], 
C(5s5pld)/[3s2pld], and H(3s)/[ls]. The atomic parameters 
are taken from the results of the X-ray structure analysis at 
room temperature. 7 ' As mentioned in ref. 16, the SOMO and 
HOMO-1 levels are well separated from the other MOs, and 
therefore, these two MOs are used to generate CAS [3,2] con- 
taining three electrons in two MOs. In order to evaluate the 
magnitudes of the inter-molecular interactions, we select two 
neighboring molecules and, in this case, put four MOs into 
the active space. Throughout the CI calculation, the MOs 
are fixed to the ones obtained from the restricted open-shell 
Hartree-Fock calculation. 

First we derive the effective model Hamiltonian describ- 
ing the isolated [TTM-TTP]+ ion. In the CAS[3,2], there 
are two independent configurations: (i) one is the ungerade 
state g 2 u* in which two electrons are on g and one electron 
is on u, and (ii) the other is the gerade state g 1 u 2 with one 
electron on g and two electrons on u. In order to evaluate 
the Coulomb interactions, we also consider the configurations 
arising from the occupations of g and u MOs by 0, 1,2 and 
4 electrons. By the CAS-CI ab initio calculations, we can ob- 
tain the full energy spectrum and the information of the wave 
functions. By collecting these data, we can construct the two- 
orbital Hubbard-type Hamiltonian that reproduces the energy 
spectrum. Namely, we consider the following one-molecule 
Hamiltonian: 

H\- mo \ —E + £g7i g + e°n u 

+ U g n Sl fn e< i + t/ u n u ,tn U 4 + U'n g n u 
1 



Jh 



S K -S U 



( c g,T c 1 L Ci a Cu <t + h - c -) 



, (2) 



where c g cr (c u cr ) is the electron annihilation operator for g 
(u) MO of the TTM-TTP molecule. The quantities e g and e° 
represent the energy levels of the g and u MOs. U s (U u ) rep- 
resents the Coulomb repulsion between two electrons on the g 
(u) MO, and U' is the inter-MO Coulomb repulsion. Jh is the 
Hund coupling including the pair-hopping term, and Eq is the 
energy constant. In contrast with the case of atomic orbitals 
under the centrosymmetric potential, there is no constraint re- 
lation among these couplings. The density operators are nor- 



Table I. Estimated parameters for the isolated TTM-TTP molecule. 
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Inter-orbital interaction 


e" = -8.63 eV 
4 = -8.21 eV 
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J H = 3.19 eV 
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Fig. 2. (Color online) Crystal structure of the (TTM-TTP)l3 projected onto 
the b—c plane (a) and onto the a—c plane (b). Definitions of transfer inte- 
grals in the MO basis (c) and in the fragment basis (d). 



mal ordered, i.e., n.„ !<T = (c^c^ a — 3/4) with v = g, u. Here, 
we note that the charges on each MO are set as 3/4 in order 
to treat the g and u MOs on the same footing. 16) The evalu- 
ated parameters are summarized in Table I. The energy level 
difference between u and g MOs is relatively small, w 0.42 
eV, which reproduces the results reported in ref. 16. The 
magnitudes of U g and U u are small compared with those of 
smaller molecules such as the TTF molecule in TTF-TCNQ, 
for which we obtained [/ttf ~ 6.2 eV. In addition, U u is 
slightly larger than U g , a reflection of the left— right bonding 
character of the u MO compared with the antibonding charac- 
ter of the g MO. The Coulomb repulsion for the TTF molecule 
was estimated from the density-functional-theory calculations 
to be sa 4.7 eV, 18) which is comparable to our result. As a 
reference, we note that the magnitude of the bare Coulomb 
repulsion of the BEDT-TTF dimer in the k-(BEDT-TTF) 2 X 
system is about 3-4 eV. 19 ' 20 ' We also note that relatively 
large Hund coupling is obtained (Jn/U u ~ 0.8), in contrast 
to those for transition-metal atoms. 

Next we evaluate the inter-molecular interactions. We fo- 
cus on two neighboring TTM-TTP molecules whose atomic 
coordinates are read from the crystal structure at room tem- 
perature, and perform CAS-CI calculations using the MOs for 
the isolated [TTM-TTP] + molecule. In order to determine this 
set of MOs, the two molecules were artificially displaced with 
respect to one another so that the intermolecular overlap was 
negligible. Hereafter, we refer to the two target molecules as 
"molecule I" and "molecule II". The molecule pairs that we 
have focused on are shown in Figs. 2(a) and 2(b). The Hamil- 
tonian for one-body terms that appear in the two-molecule 
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system is given by 



-^2-mol 



-U 



- t 



;.I.(T L u,n,CT 



^ UuC U,I,CT C U,II,CT 



<ug 



Ae g (n g j + n g , n ) 



c u,i,<x c g ,a> + h - c - 
Ae u (n u j + n u ,n) 



Ha, (3) 



where the parameters t gg and t uu are the inter-molecular 
transfer integrals for the same g and u MOs, respectively, 
while t gu and i ug are for different MOs [see Fig. 2(c)]. The 
summation over repeated spin indices is implied. In terms of 
the ab initio Hamiltonian H, these transfer integrals are ex- 
pressed as i gg = -(<p & i\H\(p S! j[), t uu = -(<p u ,i\H\<p UtS ), 
and t gu = -(ipg t i\H\ip u ^) = + (ip u j\H\ip g ^) = -t us , 
where \(p g ,i) (|<^ g ,n)) and \f u ,i) OvVn)) are the states with 
a single electron on the g MO and u MO in molecule I 
(II), respectively. The parameters Ae g and Ae u represent the 
energy-level shift due to the potential energy from the other 
molecule. The additional term is H' A = Ae gu (Cg : a c n .\. a — 

c l n cr c u,n.cr + h.c), which reflects the asymmetry between 
L/R fragments within the molecule and disappears when one 
considers the periodic crystal system. The possible inter- 
molecular Coulomb interaction terms are represented as 

-f^'-moi = ^gg^g.i^g.n + K u iujn u ,n 
+ V gn (n g ,in Ui n + n n jn g ^) 

+ I ( C g,I,(T C u,I,(T C !i,n,cr' C g:n,cr' + C g,I,cr C u,Lcr C g,n,(T' C u,n,cr') 

+ ^g( n gJ c g.n,o- c u,n,cr — c g,i,cr c u,i,o- n g,n) 



+ X u (n u ^c R .n, CT c u,n, ( 



c g,i,<T C u,i,<7 n u.n) + h - c - 



(4) 



where Vg g (V uu ) and V gu denote the inter-molecular density- 
density interactions on the same g (u) MOs and between the 
g and u MOs, respectively. The / term represents the or- 
bital exchange interaction and the X g and X u terms are the 
density-hopping interactions. As in the case of the isolated 
molecule, we can determine the parameters in eqs. (3) and 
(4) by using the ab initio energies for states with different 
symmetry, spin states, and electron numbers. The evaluated 
parameters are summarized in Table II. In contrast to the con- 
ventional extended Hiickel approach, where the overlap inte- 
grals are used in the band calculation, 2 ^ we can directly ob- 
tain the transfer integrals along different directions [ijk] in 
the present approach. We find that only the nearest-neighbor 
transfer integrals along the stacking [001] direction are large 



Table II. Estimated parameters for the inter-molecular interactions in the 
MO picture. All energies are in eV. 
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Fig. 3. (Color online) Band structure of (TTM-TTP)l3 obtained from the 
tight-binding parameters evaluated by the CAS-CI ab initio calculations 
(thick solid red lines). The dashed lines represent the energy dispersion in 
the case of t gu = 0. The GGA and its fitted results are shown by green 
dots and thin solid green lines, respectively. 



and the magnitude of ifS!? 1 ' becomes comparable to the value 
estimated by the extended Hiickel approach (« 0.26 eV). 7) 
We also note that the energy-level shifts are relatively large, 



e.g., for the [001] molecule pair, we find Ae g 001 ' = —1.81 
eV and Ael, 0011 = -1.85 eV. The magnitudes of the inter- 
molecular density-density repulsions are relatively large com- 
pared with the intra-molecular interaction (V gg /U g « 0.57 
and V uu /U u « 0.59), whereas a similar value is also obtained 
in the benchmark TTF chain in TTF-TCNQ as V/U ks 0.52. 

On the basis of the results of the above analysis, we ex- 
amine the band structure of the (TTM-TTP)l3 crystal in the 
metallic state. By neglecting correlation effects, the band 
structure obtained from the crystal version of eq. (3) is shown 
in Fig. 3. Here, we assign the energy levels for each MO as 

e g = e g + 2A4 001] = -12.25 eV and e u = e° + 2Ae l ° Q1] = 
— 11.90 eV. However, these absolute values should be re- 
fined by taking into account other molecules as well as coun- 
terion Ig\ For comparison, electronic-structure calculations 
were carried out directly for the (TTM-TTP)l3 crystal with 
the computational code QMAS (Quantum MAterials Sim- 
ulator) 22 ' based on the generalized gradient approximation 
(GGA). 23) The resulting band structure is shown in Fig. 3. De- 
spite the bandwidth being overestimated owing to the neglect 
of correlation effects, the overall band structure obtained from 
eq. (3) agrees with the GGA results, indicating the validity of 
the present approach based on the CAS-CI ab initio calcula- 
tions. The hypothetical band structure obtained by neglecting 
the SOMO-HOMO-1 mixing (i.e., we set t gu = 0) is also 
shown by the dashed lines in Fig. 3. We find that the SOMO 
band has large dispersion along the [001] direction, whereas 
the HOMO-1 band is very narrow and is located within the 
SOMO band. Furthermore, these two bands are isolated and 
thus we can conclude that this system can be regarded as a 
two-band system. 

Finally we transform the model into the fragment picture by 
using eq. (1), where the stacked TTM-TTP molecules can be 
described as an effective two-leg ladder system, as shown in 
Fig. 2(d). The transfer integral between the L and R fragments 
within the molecule is characterized by the energy difference 
between the g and u MOs, i.e., to = (e g — e u )/2 « —0.17 
eV. The Coulomb interaction within the fragment, i.e., the 
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Fig. 4. (Color online) Distance dependence of the inter-fragment Coulomb 
repulsions. The dotted line represents the bare Coulomb repulsion. The 
Coulomb repulsion within the fragment is U = 5.30 eV. 

"on-site" repulsion, is U = 5.30 eV, which becomes com- 
parable to that for the TTF molecule, Uttf ~ 6.2 eV, while 
the Coulomb and exchange interactions between the L and 
R fragments are Vb = 2.07 eV and J = 0.18 eV. Let us 
stress that Jh is not directly reduced to the exchange cou- 
pling between the L and R moieties. For the inter-molecular 
interactions, there are three kinds of transfer integrals, t\, ti, 
and £3 [shown in Fig. 2(d)], and corresponding Coulomb re- 
pulsions, Vi, V2, and V3. The correspondence relations be- 
tween the transfer integrals in the MO basis and those in the 
fragment basis are given by t\ = (— £ gg + t uu + 2t gu )/2, 
h = {tgg + i U u)/2, and t 3 = (-t gg + t uu - 2t gu )/2. The 
dominant transfer integrals (> 0.01 eV) are t = —0.17, 
4 0011 = -0.26, 4 001] = -0.17, 4 0121 = -0.05, and 4° 13] = 
0.02, where all energies are given in eV. By fitting the GGA 
band calculation, these are estimated as to = —0.12, if 101 ' = 

-0.20, if 11 = -0.14, 4 0121 = -0.03, and 4° 13] = 0.01. 
The fitted band structure is shown by the thin solid lines in 
Fig. 3. We obtain qualitatively consistent parameters, and it 
is worth noting that the inter-fragment transfer integral within 
the molecule \to\ is smaller than the inter- molecular transfer 
integral |4° 01 '|- A similar situation has been pointed out in 
single-component molecular metals. 6) 

We have also calculated the Coulomb repulsions as a func- 
tion of distance between the fragments. We define the in- 
verse of the inter-fragment distance, 1 /r, as the average in- 
verse distance between 6 sulfur atoms in each fragment. The 
magnitudes of interactions follow surprisingly well the 1/r 
Coulomb law, including the intra-molecular interaction Vb, as 
shown in Fig. 4. This result strongly supports our fragment 
decomposition picture. In conventional single-orbital com- 
pounds, the bare Coulomb interactions are known to follow 
the Coulomb law as a function of the inter-molecular dis- 
tance. 19,24) In addition, we find that the inter-fragment dis- 
tance for the [001] interactions becomes shorter than that for 
the intra-molecular one, and then the interactions v/ 001 ' and 
Vg 001 ' exceed the intra-molecular interaction Vq. This feature 
would be a key ingredient in obtaining the intra-molecular 
degree of freedom, and particularly for the ICO phenomena. 
However, we should note that the screening effects in the crys- 
tal are not taken into account. Such effects have been exam- 
ined in the single-orbital system, 18 ' 19 ' and it remains future 



work to extend the analysis to the multi-orbital systems. 

We briefly discuss the electronic states in (TTM-TTP)l3 on 
the basis of the present fragment-MO picture. From the sim- 
ple "atomic" limit analysis of our model, where the kinetic en- 
ergy is neglected, we indeed find that the lowest-energy state 
is the ICO state in which the charge is disproportionated at 
the R (L) fragment in molecule I (II) in Fig. 2(d). This ICO 
pattern is compatible with the q = (0, 0, 1/2) superstructure 
observed by the X-ray measurements. 11 ' 12, 15) On the basis of 
this finding, we infer that the non-magnetic insulating behav- 
ior observed at low temperatures' '~ 13) is attributed to the spin- 
singlet formation on the i^ 001 ' bond with two-fold periodicity 
along the stacking direction. On the other hand, the paramag- 
netic non-metallic state seen at high temperature 7-9 ' might be 
due to the charge localization on each tf 101 ' bond, since if 01 ' 
is the strongest bond in the system. A detailed analysis of pos- 
sible ordered states based on the derived effective model will 
be published elsewhere. 

In summary, we have constructed an effective Hubbard- 
type Hamiltonian and examined the electronic band structure 
for the multi-orbital molecular compound (TTM-TTP)I.3. It 
has been clarified that the present scheme combined with the 
ab initio calculations is consistent with the GGA calculation. 
We have found that, in the fragment picture, the stacked TTM- 
TTP molecules can be described by the effective two-leg lad- 
der and the Coulomb repulsions are proportional to the inverse 
of the inter-fragment distance. 
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